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In this paper we describe the off-line quality control method and 
its application in optimizing the process for forming contact windows 
in 3.5- pm complementary metal-oxide semiconductor circuits. The off- 
line quality control method is a systematic method of optimizing 
production processes and product designs. It is widely used in Japan 
to produce high-quality products at low cost. The key steps of off-line 
quality control are: (i) Identify important process factors that can be 
manipulated and their potential working levels; (ii) perform frac- 
tional factorial experiments on the process using orthogonal array 
designs; (Hi) analyze the resulting data to determine the optimum 
operating levels of the factors (both the process mean and the process 
variance are considered in this analysis; (iv) conduct an additional 
experiment to verify that the new factor levels indeed improve the 
quality control. 

I. INTRODUCTION AND SUMMARY 

This paper describes and illustrates the off-line quality control 
method, which is a systematic method of optimizing a production 
process. It also documents our efforts to optimize the process for 
forming contact windows in 3.5-jum technology complementary metal- 
oxide semiconductor (CMOS) circuits fabricated in the Murray Hill 
Integrated Circuit Design Capability Laboratory (MH ICDCL). Here, 
by optimization we mean minimizing the process variance while keep- 
ing the process mean on target. 

A typical very large scale integrated circuit (IC) chip has thousands 
of contact windows (e.g., a BELLMAC*-32 microprocessor chip has 
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250,000 windows on an approximately 1.5-cm 2 area), most of which are 
not redundant. It is critically important to produce windows of size 
very near the target dimension. (In this paper windows mean contact 
windows.) Windows that are not open or are too small result in loss of 
contact to the devices, while excessively large windows lead to shorted 
device features. The application of the off-line quality control method 
has reduced the variance of the window size by a factor of four. Also, 
it has substantially reduced the processing time required for the 
window-forming step. 

This study was inspired by Professor Genichi Taguchi's visit to the 
Quality Theory and Systems Group in the Quality Assurance Center 
at Bell Laboratories during the months of August, September, and 
October, 1980. Professor Taguchi, director of the Japanese Academy 
of Quality and a recipient of the Deming award, has developed the 
method of off-line quality control during the last three decades. It is 
used routinely by many leading Japanese industries to produce high- 
quality products at low cost. An overview of Professor Taguchi's off- 
line and on-line quality control methods is given in Taguchi, 1 and 
Kackar and Phadke. 2 This paper documents the results of the first 
application of Professor Taguchi's off-line quality control method in 
Bell Laboratories. 

The distinctive features of the off-line quality control method are 
experimental design using orthogonal arrays and the analysis of signal- 
to-noise ratios (s/n). The orthogonal array designs provide an econom- 
ical way of simultaneously studying the effects of many production 
factors on the process mean and variance. Orthogonal array designs 
are fractional factorial designs with the orthogonality property defined 
in Section IV. The s/n is a measure of the process variability. According 
to Professor Taguchi, 3 by optimizing the process with respect to the 
s/n, we ensure that the resulting optimum process conditions are 
robust or stable, meaning that they have the minimum process varia- 
tion. 

The outline of this paper is as follows: Section II gives a brief 
description of the window-forming process, which is a critical step in 
IC fabrication. The window-forming process is generally considered to 
be one of the most difficult steps in terms of reproducing and obtaining 
uniform-size windows. Nine key process factors were identified and 
their potential operating levels were determined. A description of the 
factors and their levels is given in Section III. The total number of 
possible factor-level combinations is about six thousand. 

The aim of the off-line quality control method is to determine a 
factor-level combination that gives the least variance for the window 
size while keeping the mean on target. To determine such a factor- 
level combination we performed eighteen experiments using the Lie 
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orthogonal array. The experimental setup is given in Section IV. These 
eighteen experiments correspond to eighteen factor-level combinations 
among the possible six thousand combinations. For each experiment, 
measurements were taken on the line width and the window-size 
control features. The resulting data were analyzed to determine the 
optimum factor-level combination. The measurements and the data 
analysis are presented in Sections V through IX. 

The optimum factor levels, inferred from the data analysis, were 
subsequently used in fabricating the BELLMAC-32 microprocessor, 
the BELLMAC-A microcomputer, and some other chips in the Murray 
Hill ICDCL. The experience of using these conditions is discussed in 
Section X. 

The experiment was designed and preliminary analysis of the exper- 
imental data was performed under Professor Taguchi's guidance and 
collaboration. 

II. THE WINDOW-FORMING PROCESS 

Fabrication of integrated circuits is a complex, lengthy process. 4 
Window forming is one of the more critical steps in fabricating state of 
the art CMOS integrated circuits. It comes after field and gate oxides 
are grown; polysilicon lines have been formed; and the gate, source, 
and drain areas are defined by the process of doping. Figure 1 shows 
the windows in a cross section of a wafer. A window is a hole of about 
3.5 ]Lim diameter etched through an oxide layer of about 2 /im thickness. 
The purpose of the windows is to facilitate the interconnections 
between the gates, sources, and drains. For this reason these windows 
are called contact windows. 

The process of forming windows through the oxide layers involves 
photolithography. First the P-glass surface is prepared by depositing 
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Fig. 1 — Cross section of a wafer. 
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undoped oxide on it and prebaking it. The window-forming process is 

described below. 

(i) Apply Photoresist: A wetting agent is sprayed on the wafer to 
promote adhesion of photoresist to the oxide surface. Then an appro- 
priate photoresist is applied on the wafer and the wafer is rotated at 
high speed so that the photoresist spreads uniformly. 

(ii) Bake: The wafer is baked to dry the photoresist layer. The 
thickness of the photoresist layer at this stage is about 1.3 to 1.4 /xm. 

(Hi) Expose: The photoresist-coated wafer is exposed to ultraviolet 
radiation through a mask. The windows to be printed appear as clear 
areas on the mask. In addition to the windows, which are parts of the 
desired circuits, the mask has some test patterns. Light passes through 
these areas and causes the photoresist in the window areas and the 
test pattern areas to become soluble in an appropriate solvent (devel- 
oper). The areas of the photoresist where light does not strike remain 
insoluble. 

(iv) Develop: The exposed wafer is dipped in the developer, which 
dissolves only the exposed areas. In properly printed windows, the 
exposed photoresist is removed completely and the oxide surface is 
revealed. 

(v) Plasma Etch: The wafers are placed in a high- vacuum chamber 
wherein a plasma is established. The plasma etches the exposed oxide 
areas faster than it etches the photoresist. So at the places where the 
windows are printed, windows are cut through the oxide layers down 
to the silicon surface. 

(vi) Remove Photoresist: The remaining photoresist is now re- 
moved with the help of oxygen plasma and wet chemicals. 

In the formation of the final contact windows there are additional 
steps: (vii) removal of cap-oxide, (viii) oxidation of the contact area 
to prevent diffusion of phosphorus in the subsequent step, (ix) reflow 
of the P-glass to round the window corners, (x) hydrogen annealing, 
and (xi) pre-metal wet-etching to remove any remaining oxides from 
the contact window areas. 

At the time we started this study, the target window size at step 6 
was considered to be 3.0 /mi. The final target window size (after step 
xi) was 3.5 jum. 

III. SELECTION OF FACTORS AND FACTOR LEVELS 

For the present study only the steps numbered (i) through (v) were 
chosen for optimization. Discussions with process engineers led to the 
selection of the following nine factors for controlling the window size. 
The factors are shown next to the appropriate fabrication steps. 

(i) Apply Photoresist: Photoresist viscosity (B) and spin speed 
(C). 
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(ii) Bake: Bake temperature (D) and bake time (E). 

(Hi) Expose: Mask dimension (A), aperture (F), and exposure time 
(G). 

(iv) Develop: Developing time (H). 
(v) Plasma etch: Etch time (/). 

No factor was chosen corresponding to the photoresist removal step 
because it does not affect the window size. 

The standard operating levels of the nine factors are given in Table 
I. Under these conditions, which prevailed in September 1980, the 
contact windows varied substantially in size and on many occasions 
even failed to print and open. Figure 2 shows a typical photograph of 
the programmed logic array (PLA) area of a microcomputer chip. The 
wide variation in window size and the presence of unopened windows 
is obvious from the figure. 

The principle of off-line quality control is to systematically investi- 
gate various possible levels for these factors with an aim of obtaining 
uniform-size windows. 

In the window-forming experiment a number of alternate levels were 
considered for each of the nine factors. These levels are also listed in 
Table I. Six of these factors have three levels each. Three of the factors 
have only two levels. 





Table 1 


— Test levels 






Factors Name 




Levels 


Label 




Standard 
Levels 



A Mask Dimension (fan) 2 2.5 

B Viscosity 204 206 

C Spin Speed (rpm) Low Normal High 

D Bake Temperature (°C) 90 105 

E Bake Time (min) 20 30 40 

F Aperture 12 3 

G Exposure Time 20% Over Normal 20% Under 

H Developing Time (s) 30 45 60 

I Plasma Etch Time (min) 14.5 13.2 15.8 



Dependence of spin speed on viscosity 



Spin Speed (rpm) 



Low 



Normal 



High 



Viscosity 



204 
206 



2000 
3000 



3000 
4000 



4000 
5000 



Dependence of exposure on aperture 



Exposure (PEP-Setting) 



20% Over Normal 20% Under 



Aperture 



96 
72 
40 



120 


144 


90 


108 


50 


60 
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Fig. 2— Example of nonuniform contact window sizes and an isolated, unopened 
contact window. Both are typical results obtained in August 1980, for the PLA area of 
microprocessor and microcomputer chips. The contact windows are round shaped. 



The levels of spin speed are tied to the levels of viscosity. For the 
204 photoresist viscosity the low, normal, and high spin speeds mean 
2000 rpm, 3000 rpm, and 4000 rpm, respectively. For the 206 photoresist 
viscosity the spin speed levels are 3000 rpm, 4000 rpm, and 5000 rpm. 
Likewise, the exposure setting depends on the aperture. These rela- 
tionships are also shown in Table I. 

IV. THE ORTHOGONAL ARRAY EXPERIMENT 

The full factorial experiment to explore all possible factor-level 
combinations would require 3 6 X 2 3 = 5832 experiments. Considering 
the cost of material, the time, and the availability of facilities, the full 
factorial experiment is prohibitively large. Also from statistical consid- 
erations it is unnecessary to perform the full factorial experiment 
because processes can usually be adequately characterized by a rela- 
tively few parameters. 

The fractional factorial design used for this study is given in Table 
II. It is the Lis orthogonal array design consisting of 18 experiments 
taken from Taguchi and Wu. 3 The rows of the array represent runs 
while the columns represent the factors. Here we treat BD as a joint 
factor with the levels 1, 2, and 3 representing the combinations B\D\ f 
B2D1, and B1D2, respectively. This is done so that we can study all the 

1278 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1983 



Table II — The Lie orthogonal array 



Column Number & Factor 


Experiment 


1 


2 


3 


4 


5 


6 


7 


8 


Number 


A 


BD 


C 


E 


F 


G 


H 


I 


1 




1 


1 


1 


1 


1 


1 


1 


2 




1 


2 


2 


2 


2 


2 


2 


3 




1 


3 


3 


3 


3 


3 


3 


4 




2 


1 


1 


2 


2 


3 


3 


5 




2 


2 


2 


3 


3 


1 


1 


6 




2 


3 


3 


1 


1 


2 


2 


7 




3 


1 


2 


1 


3 


2 


3 


8 




3 


2 


3 


2 


1 


3 


1 


9 




3 


3 


1 


3 


2 


1 


2 


10 


2 


1 


1 


3 


3 


2 


2 


1 


11 


2 


1 


2 


1 


1 


3 


3 


2 


12 


2 


1 


3 


2 


2 


1 


1 


3 


13 


2 


2 


1 


2 


3 


1 


3 


2 


14 


2 


2 


2 


3 


1 


2 


1 


3 


15 


2 


2 


3 


1 


2 


3 


2 


1 


16 


2 


3 


1 


3 


2 


3 


1 


2 


17 


2 


3 


2 


1 


3 


1 


2 


3 


18 


2 


3 


3 


2 


1 


2 


3 


1 



nine factors with the L\ S orthogonal array. Thus, experiment 2 would 
be run under level 1 of factors A, B, and D, and level 2 of the remaining 
factors. In terms of the actual settings, these conditions are: 2-jum mask 
dimension, 204 viscosity, 90° C bake temperature, 3000-rpm spin speed, 
bake time of 30 minutes, aperture 2, exposure PEP setting 90, 45- 
second developing time, and 13.2 minutes of plasma etch. The other 
rows are interpreted similarly. 

Here are some of the properties and considerations of this design: 
(i) This is a main-effects-only design; i.e., the response is approx- 
imated by a separable function. A function of many independent 
variables is called separable if it can be written as a sum of functions 
where each component function is a function of only one independent 
variable. 

(ii) For estimating the main effects there are two degrees of 
freedom associated with each three-level factor, one degree of freedom 
for each two-level factor, and one degree of freedom with the overall 
mean. We need at least one experiment for every degree of freedom. 
Thus, the minimum number of experiments needed is 2x6+1x3 
+ 1 = 16. Our design has 18 experiments. A single-factor-by-single- 
factor experiment would need only 16 experiments, two fewer than 18. 
But such an experiment would yield far less precise information 
compared with the orthogonal array experiment. 3,5 

(Hi) The columns of the array are pairwise orthogonal. That is, in 
every pair of columns, all combinations of levels occur and they occur 
an equal number of times. 

(iv) Consequently, the estimates of the main effects of all factors 
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as shown in Table II and their associated sums of squares are inde- 
pendent under the assumption of normality and equality of error 
variance. So the significance tests for these factors are independent. 
Though BD is treated as a joint factor, the main effects and sums of 
squares of B and D can be estimated separately under the assumption 
of no interaction. In general, these estimates would be correlated with 
each other. However, these estimates are not correlated with those for 
any of the other seven factors. 

(v) The estimates of the main effects can be used to predict the 
response for any combination of the parameter levels. A desirable 
feature of this design is that the variance of the prediction error is the 
same for all parameter-level combinations covered by the full factorial 
design. 

(vi) It is known that the main-effect-only models are liable to give 
misleading conclusions in the presence of interactions. However, in the 
beginning stages of this study the interactions are assumed to be 
negligible. If we wished to study all two-factor interactions, with no 
more than 18 experiments we would have enough degrees of freedom 
for studying only two three-level factors, or five two-level factors! That 
would mean in the present study we would have to eliminate half of 
the process factors without any experimental evidence. Alternately, if 
we wished to study all the nine process factors and their two-factor 
interactions, we would need at least 109 experiments! Orthogonal array 
designs can, of course, be used to study interactions. 3 

(vii) Optimum conditions obtained from such an experiment have 
to be verified with an additional experiment. This is done to safeguard 
us against the potential adverse effects of ignoring the interactions 
among the manipulatable factors. 

In conducting experiments of this kind, it is common for some wafers 
to get damaged or broken. Also, the wafer-to-wafer variability of 
window sizes is typically large. So we decided to run each experiment 
with two wafers. 

4.1 Analysis of variance 

Data collected from such experiments are analyzed by a method 
called analysis of variance (ANOVA). 6 The purpose of ANOVA is to 
separate the total variability of the data, which is measured by the 
sum of the squared deviations from the mean value, into contributions 
by each of the factors and the error. This is analogous to the use of 
Parseval's theorem to separate the signal strength into contributions 
by the various harmonics. 3 To see which of the factors have a signif- 
icant effect, .F-tests are performed. In performing the standard F-test 
we assume that the errors are normally distributed with equal variance 
and are independent. The results of the F-test are indicated by the 
significance level. When we say that a factor is significant at 5-percent 
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level we mean that there is 5 percent or less chance that, if we change 
the level of the factor, the response will remain the same. If the F-test 
indicates that a factor is not significant at the 5-percent level it means 
that, if we change the level of that factor, there is more than a 5- 
percent chance that the response will remain the same. 

The levels of factors which are identified as significant are then set 
to obtain the best response. The levels of the other factors can be set 
at any levels within the experimental range. We choose to leave them 
at the starting levels. 

If the assumptions of the F-test are not completely satisfied, the 
quoted significances are not accurate. However, the standard F-test is 
relatively insensitive to deviations from the assumptions used in its 
derivation. Thus, for making engineering decisions about which factor 
levels to change, the accuracy of the significance level is an adequate 
guide. In this paper we will use the standard F-test even though some 
of the assumptions are not strictly satisfied. 

V. QUALITY MEASURES 

The window size is the relevant quality measure for this experiment. 
The existing equipment does not give reproducible measurements of 
the sizes of windows in the functional circuits on a chip. This is because 
of the small size of these windows and their close proximity to one 
another. Therefore, test patterns — a line-width pattern and a window 
pattern — are provided in the upper left-hand corner of each chip. The 
following measurements were made on these test patterns to indicate 
the quality. 

(i) Line width after step (iv), called the pre-etch line width or 
photo-line width. 

(ii) Line width after step (vi), called the post-etch line width. 

(Hi) Size of the window test pattern after step (vi), called the post- 
etch window size. 

Five chips were selected from each wafer for making the above 
measurements. These chips correspond to specific locations on a 
wafer — top, bottom, left, right, and center. 

All three quality measures are considered to be good indicators of 
the size of the functional windows. However, between the geometries 
of the window-size pattern and the line-width pattern, the geometry of 
the window-size pattern is closer to the geometry of the functional 
windows. So, among the three quality measures, the post-etch window 
size may be expected to be better correlated with the size of the 
functional windows. 

VI. EXPERIMENTAL DATA 

Only thirty-four wafers were available for experimentation. So ex- 
periments 15 and 18 were arbitrarily assigned only one wafer each. One 
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Table III— Experimental data 







Line-Width Control Feature 










Photoresist — Nanoline Tool 






Experi- 






(Micrometers) 






Comments 


ment 
No. 


Top 


Center 


Bottom 


Left 


Right 


1 


2.43 


2.52 


2.63 


2.52 


2.5 




1 


2.36 


2.5 


2.62 


2.43 


2.49 




2 


2.76 


2.66 


2.74 


2.6 


2.53 




2 


2.66 


2.73 


2.95 


2.57 


2.64 




3 


2.82 


2.71 


2.78 


2.55 


2.36 




3 


2.76 


2.67 


2.9 


2.62 


2.43 




4 


2.02 


2.06 


2.21 


1.98 


2.13 




4 


1.85 


1.66 


2.07 


1.81 


1.83 




5 










— 


— 


Wafer Broke 


5 


1.87 


1.78 


2.07 


1.8 


1.83 




6 


2.51 


2.56 


2.55 


2.45 


2.53 




6 


2.68 


2.6 


2.85 


2.55 


2.56 




7 


1.99 


1.99 


2.11 


1.99 


2.0 




7 


1.96 


2.2 


2.04 


2.01 


2.03 




8 


3.15 


3.44 


3.67 


3.09 


3.06 




8 


3.27 


3.29 


3.49 


3.02 


3.19 




9 


3.0 


2.91 


3.07 


2.66 


2.74 




9 


2.73 


2.79 


3.0 


2.69 


2.7 




10 


2.69 


2.5 


2.51 


2.46 


2.4 




10 


2.75 


2.73 


2.75 


2.78 


3.03 




11 


3.2 


3.19 


3.32 


3.2 


3.15 




11 


3.07 


3.14 


3.14 


3.13 


3.12 




12 


3.21 


3.32 


3.33 


3.23 


3.10 




12 


3.48 


3.44 


3.49 


3.25 


3.38 




13 


2.6 


2.56 


2.62 


2.55 


2.56 




13 


2.53 


2.49 


2.79 


2.5 


2.56 




14 


2.18 


2.2 


2.45 


2.22 


2.32 




14 


2.33 


2.2 


2.41 


2.37 


2.38 




15 


2.45 


2.50 


2.51 


2.43 


2.43 




15 












— 


No wafer 


16 


2.67 


2.53 


2.72 


2.7 


2.6 




16 


2.76 


2.67 


2.73 


2.69 


2.6 




17 


3.31 


3.3 


3.44 


3.12 


3.14 




17 


3.12 


2.97 


3.18 


3.03 


2.95 




18 


3.46 


3.49 


3.5 


3.45 


3.57 




18 




— 


— 


— 


— 


No wafer 



of the wafers assigned to experiment 5 broke in handling. So experi- 
ments 5, 15, and 18 have only one wafer. 

The experimental data are shown in Table III. 

The data arising from such experiments can be classified as two 
types — continuous data and categorical data. Here, the pre-etch and 
the post-etch line-width data are of the continuous type. The post-etch 
window size data are mixed categorical-continuous type, because some 
windows are open while some are not. The two types of data are 
analyzed somewhat differently, as we explain the following two sec- 
tions. 

VII. ANALYSIS OF THE LINE-WIDTH DATA 

Both the pre-etch and the post-etch line widths are continuous 
variables. For each of these variables the statistics of interest are the 
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Table III — Experimental data (Continued) 







Line- Width Control Feature 










Etched— Nanoline Tool 






Experi- 






(Micrometers) 








ment 














No. 


Top 


Center 


Bottom 


Left 


Right 


Comments 


1 


2.95 


2.74 


2.85 


2.76 


2.7 




1 


3.03 


2.95 


2.75 


2.82 


2.85 




2 


3.05 


3.18 


3.2 


3.16 


3.06 




2 


3.25 


3.15 


3.09 


3.11 


3.16 




3 


3.69 


3.57 


3.78 


3.55 


3.40 




3 


3.92 


3.62 


3.71 


3.71 


3.53 




4 


2.68 


2.62 


2.9 


2.45 


2.7 




4 


2.29 


2.31 


2.77 


2.46 


2.49 




5 


— 


— 


— 


— 





Wafer Broke 


5 


1.75 


1.15 


2.07 


2.12 


1.53 




6 


3.42 


2.98 


3.22 


3.13 


3.17 




6 


3.34 


3.21 


3.23 


3.25 


3.28 




7 


2.62 


2.49 


2.53 


2.41 


2.51 




7 


2.76 


2.94 


2.68 


2.62 


2.51 




8 


4.13 


4.38 


4.41 


4.03 


4.03 




8 


4.0 


4.02 


4.18 


3.92 


3.91 




9 


3.94 


3.82 


3.84 


3.57 


3.71 




9 


3.44 


3.30 


3.41 


3.28 


3.20 




10 


3.17 


2.85 


2.84 


3.06 


2.94 




10 


3.70 


3.34 


3.45 


3.41 


3.29 




11 


4.01 


3.91 


3.92 


3.80 


3.90 




11 


3.67 


3.31 


2.86 


3.41 


3.23 




12 


4.04 


3.80 


4.08 


3.81 


3.94 




12 


4.51 


4.37 


4.45 


4.24 


4.48 




13 


3.40 


3.12 


3.11 


3.25 


3.06 




13 


3.22 


3.03 


2.89 


2.92 


2.98 




14 


3.18 


3.03 


3.4 


3.17 


3.32 




14 


3.18 


2.83 


3.17 


3.07 


3.02 




15 


2.86 


2.46 


2.3 


2.6 


2.55 




15 


— 


— 


— 


— 


— 


No wafer 


16 


2.85 


2.14 


1.22 


2.8 


3.03 




16 


3.4 


2.97 


2.96 


2.87 


2.88 




17 


4.06 


3.87 


3.90 


3.94 


3.87 




17 


4.02 


3.49 


3.51 


3.69 


3.47 




18 


4.49 


4.28 


4.34 


4.39 


4.25 




18 


— 


— 


— 


— 


— 


No wafer 



mean and the standard deviation. The objective of our data analysis 
is to determine the factor-level combination such that the standard 
deviation is minimum while keeping the mean on target. We will call 
this the optimum factor-level combination. Professor Taguchi's 
method for obtaining the optimum combination is given next. 

7.1 Single response variable 

Let us first consider the case where there is only one response 
variable. Instead of working with the mean and the standard deviation, 
it is preferable to work with the transformed variables — the mean and 
the signal-to-noise ratio (s/n). The s/n is defined as 

\ Standard Deviation/ 
= — logio(coefficient of variation). 
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Table III— Experimental data (Continued) 







Window-Control Feature 








Etched 


— Vickers Tool 




Experi- 




(Micrometers) 






ment 
No. 


Top 


Center 


Bottom 


Left 


Right Comments 


1 


WNO* 


WNO 


WNO 


WNO 


WNO 


1 


WNO 


WNO 


WNO 


WNO 


WNO 


2 


2.32 


2.23 


2.30 


2.56 


2.51 


2 


2.22 


2.33 


2.34 


2.15 


2.35 


3 


2.98 


3.14 


3.02 


2.89 


3.16 


3 


3.15 


3.08 


2.78 


WNO 


2.86 


4 


WNO 


WNO 


WNO 


WNO 


WNO 


4 


WNO 


WNO 


WNO 


WNO 


WNO 


5 










— 


— Wafer Broke 


5 


WNO 


WNO 


WNO 


WNO 


WNO 


6 


2.45 


2.19 


2.14 


2.32 


2.12 


6 


WNO 


WNO 


WNO 


WNO 


WNO 


7 


WNO 


WNO 


WNO 


WNO 


WNO 


7 


WNO 


WNO 


WNO 


WNO 


WNO 


8 


WNO 


WNO 


WNO 


WNO 


WNO 


8 


2.89 


2.97 


3.13 


3.25 


3.19 


9 


3.16 


2.91 


3.12 


3.18 


3.11 


9 


2.43 


2.35 


2.14 


2.40 


2.28 


10 


2.0 


1.75 


1.97 


1.91 


1.72 


10 


WNO 


2.7 


WNO 


2.61 


2.73 


11 


2.76 


3.09 


3.22 


3.05 


3.04 


11 


3.12 


3.21 


WNO 


2.71 


2.27 


12 


3.24 


3.08 


WNO 


2.89 


2.72 


12 


3.5 


3.71 


3.52 


3.53 


3.71 


13 


2.54 


2.63 


2.88 


2.31 


2.71 


13 


WNO 


WNO 


WNO 


WNO 


WNO 


14 


WNO 


1.74 


2.24 


2.07 


2.38 


14 


WNO 


WNO 


WNO 


WNO 


WNO 


15 


WNO 


WNO 


WNO 


WNO 


WNO 


15 










— 


— No wafer 


16 


WNO 


WNO 


WNO 


WNO 


WNO 


16 


WNO 


WNO 


WNO 


WNO 


WNO 


17 


3.09 


2.91 


3.06 


3.09 


3.29 


17 


3.39 


2.5 


2.57 


2.62 


2.35 


18 


3.39 


3.34 


3.45 


3.44 


3.33 


18 






— 


— 


— No wafer 



* WNO— Window not open. 

In terms of the transformed variables, the optimization problem is 
to determine the optimum factor levels such that the s/n is maximum 
while keeping the mean on target. This problem can be solved in two 
stages: 

(i) Determine which factors have a significant effect on the s/n. 
This is done through the analysis of variance (ANOVA) of the s/n. 
These factors are called the control factors, implying that they control 
the process variability. For each control factor we choose the level 
with the highest s/n as the optimum level. Thus the overall s/n is 
maximized. 

(«') Select a factor that has the smallest effect on the s/n among all 
factors that have a significant effect on the mean. Such a factor is 
called a signal factor. Ideally, the signal factor should have no effect 
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on the s/n. Choose the levels of the remaining factors (factors that are 
neither control factors nor signal factors) to be the nominal levels prior 
to the optimization experiment. Then set the level of the signal factor 
so that the mean response is on target. 

In practice, the following two aspects should also be considered in 
selecting the signal factor: (i) If possible, the relationship between the 
mean response and the levels of the signal factor should be linear, and 
(u) It should be convenient to change the signal factor during produc- 
tion. These aspects are important from the on-line quality control 
considerations. The signal factor can be used during manufacturing to 
adjust the mean response. 1 " 3 

Why do we work in terms of the s/n ratio rather than the standard 
deviation? Frequently, as the mean decreases, the standard deviation 
also decreases and vice versa. In such cases, if we work in terms of the 
standard deviation, the optimization cannot be done in two steps; i.e., 
we cannot minimize the standard deviation first and then bring the 
mean on target. 

Through many applications, Professor Taguchi has empirically 
found that the two-stage optimization procedure involving the s/n 
indeed gives the parameter-level combination where the standard 
deviation is minimum, while keeping the mean on target. This implies 
that the engineering systems behave in such a way that the manipu- 
latable production factors can be divided into three categories: 

(i) Control factors, which affect process variability as measured by 
the s/n 

(ii) Signal factors, which do not influence (or have negligible effect 
on) the s/n but have a significant effect on the mean 

(Hi) Factors that do not affect the s/n or the process mean. 

The two-stage procedure also has an advantage over a procedure 
that directly minimizes the mean square error from the target mean 
value. In practice, the target mean value may change during the 
process development. The advantage of the two-stage procedure is 
that for any target mean value (of course, within reasonable bounds) 
the new optimum factor-level combination is obtained by suitably 
adjusting the level of only the signal factor. This is so because in step 
(i) of the algorithm the coefficient of variation is minimized for every 
mean target value. 

7.2 Multiple response variables 

Now let us consider the case where there are two or more response 
variables. In such cases, engineering judgment may have to be used to 
resolve the conflict if different response variables suggest different 
levels for any one factor. The modified two-stage procedure is as 
follows: 
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(i) Separately determine control factors and their optimum levels 
corresponding to each response variable. If there is a conflict between 
the optimum levels suggested by the different response variables, use 
engineering judgment to resolve the conflict. 

(ii) Select a factor that has the smallest effect (preferably no effect) 
on the signal-to-noise ratios for all the response variables but has a 
significant effect on the mean levels. This is the signal factor. Set the 
levels of the remaining factors, which affect neither the mean nor the 
s/n, at the nominal levels prior to the optimization experiment. Then 
set the level of the signal factor so that the mean responses are on 
target. Once again engineering judgment may have to be used to 
resolve any conflicts that arise. 

The selection of the control factors, signal factor, and their optimum 
levels for the present application will be discussed in Section IX. The 
remaining portions of Sections VII and VIII contain the data analysis 
that forms the basis for selecting the optimum factor levels. 

7.3 Pre-etch line width 

Mean, standard deviation, and s/n were calculated for each of the 
eighteen experiments. For those experiments with two wafers, ten data 
points were used in these calculations. When there was only one wafer, 
five data points were used. These results are shown in Table IV. The 
presence of unequal sample sizes has been ignored in the subsequent 
analysis. Let Xi and tj, denote the mean and the s/n for the ith 
experiment. 

Table IV— Pre-etch line-width data 







Standard 








Deviation 






Mean 


of 






Line Width, 


Line Width, 




Experiment 


* 


8 


s/n 


Number 


(/un) 


(lim) 


tj = log(x/s) 


1 


2.500 


0.0827 


1.4803 


2 


2.684 


0.1196 


1.3512 


3 


2.660 


0.1722 


1.1889 


4 


1.962 


0.1696 


1.0632 


5 


1.870 


0.1168 


1.2043 


6 


2.584 


0.1106 


1.3686 


7 


2.032 


0.0718 


1.4520 


8 


3.267 


0.2101 


1.1917 


9 


2.829 


0.1516 


1.2709 


10 


2.660 


0.1912 


1.1434 


11 


3.166 


0.0674 


1.6721 


12 


3.323 


0.1274 


1.4165 


13 


2.576 


0.0850 


1.4815 


14 


2.308 


0.0964 


1.3788 


15 


2.464 


0.0385 


1.8065 


16 


2.667 


0.0706 


1.5775 


17 


3.156 


0.1569 


1.3036 


18 


3.494 


0.0473 


1.8692 
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By computing a single mean x t and a single variance sf (needed for 
computing tj,-) from the two wafers of each experiment i, we pool 
together the between wafer and the within wafer variance. That is, 

E(sf) = (Between wafer variance for experiment i) x % 

+ (Within wafer variance for experiment i). 

Thus, when we maximize % we minimize the sum of the between- wafer 
and the within-wafer variances of the line width, which is the response 
of interest to us. There can be situations when one wants to separately 
estimate the effects of the factor levels on the between-wafer and 
within-wafer variances. In those cases, one would compute the s/n and 
the mean line width for each individual wafer. 

In the analysis of the pre-etch and the post-etch line widths, we 
compute the x, and the sf for each experiment by pooling the data 
from both wafers used in that experiment. A relative measure of the 
between-wafer and within-wafer variance is obtained in Section VIII, 
while the post-etch window-size data is being analyzed. 

7.3.1 Analysis of s/n 

The estimates of the average s/n for all factor levels are given in 
Table V. The average for the first level of factor A is the average of 
the nine experiments (experiments 1 through 9), which were conducted 
with level 1 of the factor A. Likewise, the average for the second level 
of factor A is the mean of experiments 10 through 18, which were 
conducted with level 2 of the factor A. Let us denote these average 
effects of Ai and A 2 by ntA 1 and tha 2 , respectively. Here m.A l = 1.2857 
and m,A 2 = 1.5166. The other entries of Table V were calculated 
similarly. 

The average signal-to-noise ratios for every level of the eight factors 
are graphically shown in Fig. 3. Qualitatively speaking, the mask 
dimension and the aperture cause a large variation in the s/n. The 



Table V — Pre-etch line width for average s/n 







Average s/n 




Factor 


Level 1 


Level 2 


Level 3 


A Mask Dimension 


1.2857 


1.5166 




BD Viscosity Bake Temperature 


(BiDi) 1.3754 


(fi 2 D,) 1.3838 


(BiD 2 ) 1.4442 


B Viscosity 


1.4098 


1.3838 




D Bake Temperature 


1.3796 


1.4442 




C Spin Speed 


1.3663 


1.3503 


1.4868 


E Bake Time 


1.4328 


1.4625 


1.3082 


F Aperture 


1.5368 


1.4011 


1.2654 


G Exposure Time 


1.3737 


1.3461 


1.4836 


H Developing Time 


1.3881 


1.4042 


1.4111 



Overall average s/n = 1.4011. 
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developing time and the viscosity cause a small change in the s/n. The 
effect of the other factors is in between. 

For three-level factors, Fig. 3 can also be used to judge the linearity 
of the effect of the factors. If the difference between levels 1 and 2, and 
levels 2 and 3 is equal and these levels appear in proper order (1, 2, 3, 
or 3, 2, 1), then the effect of that factor is linear. If either the differences 
are unequal or the order is mixed up, then the effect is not linear. For 
example, the aperture has approximately linear response while the 
bake time has a nonlinear response. 

We shall perform a formal analysis of variance (ANOVA) to identify 
statistically significant factors. The analysis of variance of general 
linear models is widely known in literature, see e.g., Searle 7 and Hicks. 6 
Simple ANOVA methods for orthogonal array experiments are de- 
scribed in Taguchi and Wu. 3 The linear model used in analyzing this 
data is: 

yi = n + Xi + et, (1) 

where 

i = 1, • • • , 18 is the experiment number. 

li is the overall mean. 

Xi is the fixed effect of the factor-level combination used in exper- 
iment i. Here we consider only the main effect for each of the 
factors. Thus it represents the sum of the effects of the eight 
factors. 

e, is the random error for experiment i. 

yi is the s/n for experiment i. 

To clarify the meaning of the term Xi, let us consider experiment 1, 
which was run at level 1 of each of the eight factors A through H. Note 
that the factor I is irrelevant for studying the pre-etch line width. So 
Xi is the sum of the main effects associated with the first level of each 
of the factors A through H. 

The sum of squares and the mean squares for the eight factors are 
tabulated in Table Via. The computations are illustrated in Appendix 
A. 

The expected mean squares are also shown in Table Via. See Refs. 
6 and 7 for the computation of expected mean squares, which are used 
in forming appropriate F- tests. The error variance, i.e., variance of e,, 
is denoted by a 2 . The variability due to the factors A through H is 
denoted by <$> with an appropriate subscript. 

In Table Via we see that the mean sum of squares for factors BD, 
C, G, and H are smaller than the mean error sum of squares. So a new 
ANOVA table, Table VIb, was formed by pooling the sum of squares 
of these factors with the error sum of squares. The linear model 
underlying the ANOVA Table VIb is the same as Eq. (1), except that 
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now Xi stands for the sum of the main effects of only A, E, and F. The 
F ratios, computed by dividing the factor mean square by the error 
mean square, are also shown in Table VIb. Factors A and F are 
significant using F-table values for the 5-percent significance level. So 
the mask dimension and the aperture are the control factors. 

In performing the analysis of variance, we have tacitly assumed that 
the response for each experiment, here the s/n, has a normal distri- 
bution with constant variance. We are presently investigating the 
distributional properties of the s/n and their impact on the analysis of 
variance. In this paper we treat the significance levels as approximate. 

The engineering significance of a statistically significant factor can 
be measured in terms of the percent contribution, a measure intro- 
duced by Taguchi. 3 The percent contribution is equal to the percent of 
the total sum of squares explained by that factor after an appropriate 
estimate of the error sum of squares has been removed from it. The 
larger the percent contribution, the more can be expected to be 
achieved by changing the level of that factor. Computation of the 
percent contribution is illustrated in Appendix B, and the results are 
shown in Table VIb. 

From Table VIb we see that both the factors A (mask dimension) 
and F (aperture) contribute in excess of 20 percent each to the total 



Table VI — Pre-etch line width 



(a) ANOVA for s/n 


Source 




Degrees 
of Free- 
dom 


Sum of 
Squares 


Mean 
Square 


Expected 
Mean 
Square 


A Mask Dimension 1 
BD Viscosity Bake Temperature 2 
C Spin Speed 2 
E Bake Time 2 
F Aperture 2 
G Exposure Time 2 
H Developing Time 2 
Error 4 

Total 17 


0.2399 
0.0169 
0.0668 
0.0804 
0.2210 
0.0634 
0.0017 
0.1522 

0.8423 


0.2399 
0.0085 
0.0334 
0.0402 
0.1105 
0.0317 
0.0009 
0.0381 


C? + <t>A 
O + <j>BD 

a 2 + <t>c 

t^ + ^E 

a* + <t> F 
a 2 + 4> c 

2 + <f>H 
CT 2 


(b) Pooled ANOVA for s/n 


Source 


Degrees 
of Free- 
dom 


Sum of 
Squares 


Mean 
Square 


F 


Percent 
Contribution 


A Mask Dimension 
E Bake Time 
F Aperture 
Error 

Total 


1 
2 

2 
12 

17 


0.2399 
0.0804 
0.2210 
0.3010 

0.8423 


0.2399 
0.0402 
0.1105 
0.0251 


9.56* 

1.60 

4.40* 


25.5 

3.6 

20.3 

50.6 

100.00 



F,.,2(0.95) = 4.75. 
F 2 .i 2 (0.95) = 3.89. 
* Factors significant at 95-percent confidence level. 
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sum of squares. So the factors A and F are not only statistically 
significant, they have a sizable influence on the s/n. These results are 
consistent with Fig. 3. They will be used in Section IX for selecting the 
control factors. 

7.3.2 Analysis of the means 

Now we analyze the mean pre-etch line widths, Xi values, to find a 
signal factor. 

The estimates of the mean line widths for all factor levels are given 
in Table VII. These estimates are graphically shown in Fig. 4. It is 
apparent that the levels of viscosity, mask dimension, and spin speed 
cause a relatively large change in the mean line width. Developing 
time and aperture have a small effect on the line width. The remaining 
two factors have an intermediate effect. 

The linear model used to analyze this data is the same as eq. (1), 
except that now yi stands for the mean pre-etch line width rather than 
the s/n. 

The original and the pooled ANOVA tables for the mean pre-etch 
line width are given in Tables Villa and b, respectively. Because the 
design is not orthogonal with respect to the factors B and D, we need 
a special method, described in Appendix C, to separate Sbd into Sb 
and Sd. 

It is clear from Table VHIb that the mask dimension (A), viscosity 
(B), and spin speed (C) have a statistically significant effect on the 
mean pre-etch line width. Also, these factors together contribute more 
than 70 percent to the total sum of squares. These results will be used 
in Section IX for selecting the signal factor. 

7.4 Post-etch line width 

The analysis of the post-etch line-width data is similar to the 
analysis of the pre-etch line-width data. The mean, the standard 
deviation, and the s/n for each experiment are shown in Table IX. 

Table VII — Pre-etch line width for the mean line width 



Mean Line Width (/im) 


Factor 


Level 1 


Level 2 


Level 3 


A Mask Dimension 


2.39 


2.87 




BD Viscosity Bake Temperature 


(B,D,) 2.83 


(BiDx) 2.31 


(5,D 2 ) 2.74 


B Viscosity 


2.79 


2.31 




D Bake Temperature 


2.57 


2.74 




C Spin Speed 


2.40 


2.59 


2.89 


E Bake Time 


2.68 


2.68 


2.53 


F Aperture 


2.68 


2.56 


2.64 


G Exposure Time 


2.74 


2.66 


2.49 


H Developing Time 


2.60 


2.60 


2.69 



Overall mean line width = 2.63 |im. 
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Table VIM— Pre-etch line width 



(a) ANOVA for mean line width 



Source 



Degrees 
of Free- 
dom 



Sum of 
Squares 



Mean 
Square 



Expected 
Mean 
Square 



A Mask Dimension 

BD Viscosity Bake Temperature 

C Spin Speed 

E Bake Time 

F Aperture 

G Exposure Time 

H Developing Time 

Error 



1 


1.05 


1.050 


2 


0.95 


0.475 


2 


0.73 


0.365 


2 


0.10 


0.050 


2 


0.05 


0.025 


2 


0.19 


0.095 


2 


0.04 


0.020 


4 


0.26 


0.065 



A Mask Dimension 


1 


1.05 


1.050 


19.81* 


B Viscosity 


1 


0.83 


0.834 


15.74* 


C Spin Speed 


2 


0.73 


0.365 


6.89* 


G Exposure Time 


2 


0.19 


0.095 


1.79 


Error 


11 


0.58 


0.053 





Total 



17 



3.37 



£ + ** 

» + 't'BD 

<r + <j> E 
a 2 + fa 
a 



Total 




17 3.37 




(b) Pooled ANOVA for mean line width 




Source 


Degrees 

of Free- Sum of Mean 
dom Squares Square F 


Contribution 
Percent 



29.6 
22.6 
18.5 
2.5 
26.8 

100.0 



Fi.ii (0.95) = 4.84. 
F 2 .n(0.95) = 3.98. 
* Factors significant at 95-percent confidence level. 



Table IX— Post-etch line-width data 







Standard 






Mean 


Deviation of 






Line Width, 


Line Width 




Experiment 


X 


s 


s/n 


Number 


{fim) 


(/tm) 


t; = \og{x/s) 


1 


2.84 


0.11 


1.42 


2 


3.14 


0.063 


1.70 


3 


3.65 


0.15 


1.40 


4 


2.57 


0.20 


1.11 


5 


1.72 


0.40 


0.63 


6 


3.12 


0.27 


1.07 


7 


2.62 


0.19 


1.14 


8 


4.10 


0.18 


1.37 


9 


3.55 


0.26 


1.13 


10 


3.31 


0.35 


0.98 


11 


3.60 


0.38 


0.98 


12 


4.17 


0.27 


1.18 


13 


3.10 


0.16 


1.29 


14 


3.14 


0.16 


1.29 


15 


2.55 


0.21 


1.09 


16 


2.81 


0.37 


0.88 


17 


3.78 


0.22 


1.23 


18 


4.34 


0.078 


1.75 
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The average s/n and the mean line width for each factor level are 
shown in Tables Xa and b, respectively. 

The linear model (1) was again used to analyze the post-etch line- 
width data. The ANOVA for the signal-to-noise ratios, Table XIa, 
indicates that none of the nine process factors has a significant effect 
(approximately 5-percent level) on the s/n for the post-etch line width. 
The pooled ANOVA for the mean post-etch line widths is shown in 
Table Xlb. It is obvious from the table that the viscosity, exposure, 
spin speed, mask dimension, and developing time have significant 
effects (5-percent level) on the mean line width. The contribution of 
these factors to the total sum of squares exceeds 90 percent. The mean 
line width for each factor level is shown graphically in Fig. 5. 

VIII. ANALYSIS OF POST-ETCH WINDOW-SIZE DATA 

Some windows are printed and open while the others are not. Thus 
the window-size data are mixed categorical-continuous in nature. Anal- 
ysis of such data is done by converting all the data to the categorical 
type and then using the 'accumulation analysis' method, which is 



Table X— Post-etch line width 





(a) Average signal-to-noise 


ratios 












Average s/n 




Factor 


Level 1 


Level 2 


Level 3 


A Mask Dimension 




1.22 




1.19 




BD Viscosity Bake Temperature 


(BiDi) 1.28 


(BzA) 1.08 


(BiZh) 1.25 


B Viscosity 




1.27 




1.08 




D Bake Temperature 




1.18 




1.25 




C Spin Speed 




1.14 




1.20 


1.27 


E Bake Time 




1.16 




1.28 


1.17 


F Aperture 




1.28 




1.22 


1.11 


G Exposure Time 




1.26 




1.33 


1.02 


H Developing Time 




1.09 




1.20 


1.32 


I Etch Time 




1.21 




1.18 


1.23 



Overall average s/n = 1.205 



(b) Mean line width 



Mean Line Width (/xm) 



Factor 


Level 1 


Level 2 


Level 3 


A Mask Dimension 


3.03 


3.42 




BD Viscosity Bake Temperature 


(£iA) 3.45 


(£2Di ) 2.70 


(£,Z) 2 ) 3.53 


B Viscosity 


3.49 


2.70 




D Bake Temperature 


3.08 


3.53 




C Spin Speed 


2.88 


3.25 


3.56 


E Bake Time 


3.15 


3.18 


3.35 


F Aperture 


3.28 


3.22 


3.18 


G Exposure Time 


3.52 


3.34 


2.83 


H Developing Time 


3.04 


3.09 


3.56 


I Etch Time 


3.14 


3.22 


3.32 



Overall mean line width = 3.23 fan. 
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described by Taguchi in Refs. 3 and 8. Factors that are found significant 
in this analysis are control factors. 

The window sizes were divided into the following five categories: 



Category 



Description (micrometers) 



I 


Window not open or not printed 


II 


(0, 2.25) 


III 


[2.25, 2.75) 


IV 


[2.75, 3.25] 


V 


(3.25, oo) 



Note that these categories are ordered with respect to window size. 
The target window size at the end of step (vi) was 3 /xm. Thus category 
IV is the most desired category, while category I is the least desired 
category. Table XII summarizes the data for each of the experiments 
by categories. To simplify our analysis, we shall presume that a missing 
wafer has the same readings as the observed wafer for that experiment. 
This is reflected in Table XII, where we show the combined readings 
for the two wafers of each experiment. 

Table XI— Post-etch line width 

a) ANOVA for s/n 



Source 



Degrees 
of Free- 
dom 



Sum of 
Squares 



Mean 

Square 



A Mask Dimension 


1 


0.005 


0.005 


0.02 


B Viscosity 


1 


0.134 


0.134 


0.60 


D Bake Temperature 


1 


0.003 


0.003 


0.01 


C Spin Speed 


2 


0.053 


0.027 


0.12 


E Bake Time 


2 


0.057 


0.028 


0.13 


F Aperture 


2 


0.085 


0.043 


0.19 


G Exposure Time 


2 


0.312 


0.156 


0.70 


H Developing Time 


2 


0.156 


0.078 


0.35 


I Etch Time 


2 


0.008 


0.004 


0.02 


Error 


2 


0.444 


0.222 





Total 


17 1.257 




b) Pooled ANOVA for mean line width 


Source 


Degrees 

of Free- Sum of Mean 
dom Squares Square F 


Percent 
Contri- 
bution 



A Mask Dimension 


1 


0.677 


0.677 


16.92* 


8.5 


B Viscosity 


1 


2.512 


2.512 


63.51* 


32.9 


C Spin Speed 


2 


1.424 


0.712 


17.80* 


17.9 


G Exposure Time 


2 


1.558 


0.779 


19.48* 


19.6 


H Developing Time 


2 


0.997 


0.499 


12.48* 


12.2 


Error 


9 


0.356 


0.040 




8.9 


Total 


17 


7.524 






100.0 



F 1>9 (0.95) = 5.12. 
F 2 . 9 (0.95) = 4.26. 
* Factors significant at 95-percent confidence level. 
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Table XII — Post-etch window-size data — frequencies by experiment 





Freq 


[uency Distribution 


Freq 


[uency Distribution 


Combined Frequency for 


Experi- 




for Wafer 1 






for Wafer 2 






the Two Wafers 




ment 
































No. 


I 


n 


III 


IV 


V 


I 


n 


III 


IV 


V 


I 


n 


m 


IV 


V 


1 


5 














5 














10 














2 





1 





2 


2 





2 


3 











3 


3 


2 


2 


3 











4 





1 








3 





1 








9 





4 


5 














5 














10 














5 


* 


• 


* 


* 


* 


5 














10 














6 





3 


2 








5 














5 


3 


2 








7 


5 














5 














10 














8 


5 























5 





5 








5 





9 











5 








1 


4 











1 


4 


5 





10 





5 











2 





3 








2 


5 


3 








11 











5 





1 


1 


2 


1 





1 


1 


2 


6 





12 


1 





1 


3 

















4 


1 





1 


3 


5 


13 








3 


2 





5 














5 





3 


2 





14 


1 


3 


1 








5 














6 


3 


1 








15 


5 














* 


* 


» 


* 


« 


10 














16 


5 














5 














10 














17 











3 


2 








4 





1 








4 


3 


3 


18 














5 


» 


» 


* 


* 


* 














10 



* Implies data missing. 

Table XIII gives the frequency distribution corresponding to each 
level of each factor. To obtain the frequency distribution for a specific 
level of a specific factor, we summed the frequencies of all the experi- 
ments that were conducted with that particular level of that particular 
factor. For example, the frequency distribution for the first level of 
factor C (low spin speed) was obtained by summing the frequency 
distributions of experiments with serial numbers 1, 4, 7, 10, 13, and 16. 
These six experiments were conducted with level 1 of factor C. 

The frequency distributions of Table XIII are graphically displayed 
by star plots in Fig. 6. From this figure and the table it is apparent 
that a change in the level of viscosity, spin speed, or mask dimension 
causes a noticeable change in the frequency distribution. A change in 
the level of etch time, bake time, or bake temperature seems to have 
only a small effect on the frequency distribution. The effects of the 
other factors are intermediate. 

We now determine which factors have a significant effect on the 
frequency distribution of the window sizes. The standard chi-square 
test for multinomial distributions is not appropriate here because the 
categories are ordered. The accumulation analysis method has an 
intuitive appeal and has been empirically found by Professor Taguchi 
to be effective in analyzing ordered categorical data. The method 
consists of the following three steps: 

(i) Compute the cumulative frequencies. Table XIII shows the 
cumulative frequencies for all factor levels. The cumulative categories 
are denoted with parentheses. Thus (III) means sum of categories I, 

IC QUALITY CONTROL 1297 



Table XIII— Post-etch window-size data— frequencies by factor level 







Frequencies 






Cumulative Frequencies 




Factor Levels 


I 


n 


m 


IV 


V 


(I) 


(ID 


(HI) 


(IV) 


(V) 


Mask Dimension 






















A, 


51 


7 


9 


21 


2 


51 


58 


67 


88 


90 


A 2 


35 


9 


14 


14 


18 


35 


44 


58 


72 


90 


Viscosity, Bake 






















Temperature 






















BtDx 


15 


9 


9 


20 


7 


15 


24 


33 


53 


60 


B2D1 


46 


6 


6 


2 





46 


52 


58 


60 


60 


BiDi 


25 


1 


8 


13 


13 


25 


26 


34 


47 


60 


Spin Speed 






















G 


47 


5 


6 


2 





47 


52 


58 


60 


60 


c 2 


22 


7 


10 


16 


5 


22 


29 


39 


55 


60 


C 3 


17 


4 


7 


17 


15 


17 


21 


28 


45 


60 


Bake Time 






















£1 


31 


2 


10 


14 


3 


31 


33 


43 


57 


60 


E 2 


26 


3 


7 


7 


17 


26 


29 


36 


43 


60 


hiz 


29 


11 


6 


14 





29 


40 


46 


60 


60 


Aperture 






















F, 


32 


7 


5 


6 


10 


32 


39 


44 


50 


60 


F 2 


36 


3 


4 


10 


7 


36 


39 


43 


53 


60 


F a 


18 


6 


14 


19 


3 


18 


24 


38 


57 


60 


Exposure Time 






















26 


3 


10 


13 


8 


26 


29 


39 


52 


60 


G 2 


18 


12 


11 


7 


12 


18 


30 


41 


48 


60 


G 3 


42 


1 


2 


15 





42 


43 


45 


60 


60 


Developing Time 






















H x 


37 


4 


6 


8 


5 


37 


41 


47 


55 


60 


H 2 


27 


11 


12 


5 


5 


27 


38 


50 


55 


60 


H 3 


22 


1 


5 


22 


10 


22 


23 


28 


50 


60 


Etch Time 






















/1 


37 


5 


3 


5 


10 


37 


42 


45 


50 


60 


h 


21 


8 


14 


15 


2 


21 


29 


43 


58 


60 


h 


28 


3 


6 


15 


8 


28 


31 


37 


52 


60 


Totals 


86 


16 


23 


35 


20 


86 


102 


125 


160 


180 



II, and III. Note that the cumulative category (V) is the same as the 
total number of window readings for the particular factor level. 

(ii) Perform "binary data" ANOVA 7 on each cumulative category 
except the last category, viz. (V). Note that a certain approximation is 
involved in the significance levels suggested by this ANOVA because 
the observations are not normally distributed. 

{Hi) Assign weights to each cumulative category. These weights are 
inversely proportional to the Bernoulli trial variance. Let cum c be the 
total number of windows in the cumulative category, c, as given in the 
bottom row of Table XIII. Then the weight for that category is: 

w.- 1 18 ° 2 



cum c ( cum c \ cum c (180 - cum c ) 
x 1 1 — 



180 V 180 

These weights are shown in Appendix D for each category. 
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Then for each factor and for each error term the accumulated sum 
of squares is taken to be equal to the weighted sum of the sum of 
squares for all cumulative categories. 

The intuitive appeal for accumulation analysis is that by taking 
cumulative frequencies we preserve the order of the categories. By 
giving weights inversely proportional to the sampling errors in each 
cumulative category, we make the procedure more sensitive to a 
change in the variance. The difficulty is that the frequencies of the 
cumulative categories are correlated. So the true level of significance 
of the F-test may be somewhat different from that indicated by the F 
table. More work is needed to understand the statistical properties of 
the accumulation analysis. 

Table XIV gives the final ANOVA with accumulated sum of squares. 
The computations are illustrated in Appendix D. For each cumulative 
category, the following nested, mixed linear model was used in per- 
forming the ANOVA: 

yijk = n + Xi + euj + e 2 ijk, (2) 

where 

i = 1, • • • , 18 stands for the experiment. 

j = 1, 2 stands for wafer within the experiment. 

k = 1, • • • , 5 stands for replicate or position within wafer within 

the experiment, 
/i is the overall mean. 

Xi is the fixed effect of the factor-level combination used in 

experiment i. Here we consider only the main effect for each 

of the factors. See the discussion of model (1) in Section 7.3 

for more details of the interpretation of Xi. 

euj is the random effect for wafer j within experiment i. 

e%jjK is the random error for replicate k within wafer j within 

experiment i. 
yijk is the observation for replicate k in wafer j in experiment i. 
yijk takes a value 1 if the window size belongs to the particular 
category. Otherwise, the value is zero. 
The expected mean squares for this ANOVA model are also shown 
in Table XIV. The variances of e\ and e 2 are denoted by a\ and a\, 
respectively. The effects of the factors A through I are denoted by <f> 
with an appropriate subscript. The effect of lack of fit is denoted by 
<f>s. We assume that the random variables euj and e2ijk are independent 
for all values of i,j and k. The degrees of freedom shown in Table XIV 
have been adjusted for the fact that three experiments have only one 
wafer each. 

For testing the significance of the effect of error between wafers 
within experiments, the relevant denominator sum of squares is the 
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Table XIV — Post-etch window size 



ANOVA for accumulation analysis 




Degrees 










of Free 


i- Sum of 


Mean 




(Expected Mean 


Source 


dom 


Squares 


Square 


F 


Square) + W 


A Mask Dimension 


4 


26.64 


6.66 


2.67* 


&2 + 5(Ti + <j>A 


BD Viscosity-Bake 


8 


112.31 


14.04 


5.64* 


02 + 501 + <fiBD 


Temperature 












C Spin Speed 


8 


125.52 


15.69 


6.30* 


O-l + 5ffi 2 + $C 


E Bake Time 


8 


36.96 


4.62 


1.86 


a-i + 5ui 4- <f>E 


F Aperture 


8 


27.88 


3.49 


1.40 


Ol + 5ffi + <f>F 


G Exposure Time 


8 


42.28 


5.29 


2.12* 


<J2 + 5(Ti + <f>G 


H Developing Time 


8 


45.57 


5.70 


2.29* 


02 + 5di + <j>H 


I Etch Time 


8 


23.80 


2.98 


1.20 


Oi + 5ffi + <f>l 


Lack of Fit 


8 


17.25 


2.16 


0.87 


02 + 5ai + <f>i 


Error Between Wafers 


60 


149.33 


2.49 


11.69* 


02 2 + 5ai 2 


Within Experiment 












Error Between Repli- 


528 


112.45 


0.21 




Of 


cates Within Wafers 












Within Experiment 












Total 


656 


720.00 








W- (Wo) + W m + Wum + W (lV) )/4. 








F4,6o(0.95) = 2.53, F 8 .6o 


(0.95) = 


2.10, F6o.528(0.95) = 1.32. 
(b) Separation of Sbd 










Degrees 












of Free- 


Sum of 


Mean 




Source 




dom 


Squares 


Square 


F 


B Viscosity 




4 


87.38 


21.85 


8.78* 


D Bake Temperature 


4 


6.55 


1.64 


0.66 



* Factors significant at 95-percent confidence level. 



estimate of ah The corresponding F value is 11.69, which is significant 
far beyond the nominal 5-percent level. To test for the lack of fit of the 
main-effects-only model, the appropriate denominator is the estimate 
ai + 5a i. The corresponding F ratio is 0.87. This indicates that the 
main-effects-only model adequately describes the observed data rela- 
tive to the random errors between wafers. For testing the significance 
of the process factors, the denominator mean square is again the 
estimate of a\ + 5a?. We see that the mask dimension, viscosity, spin 
speed, exposure time, and developing time have a significant effect 
(approximately 5-percent level) on the window size. The effects of the 
other factors are not significant. 

IX. SELECTION OF OPTIMUM FACTOR LEVELS 

The following table summarizes the significant results of the analyses 
performed in Sections VII and VIII. In each category, the factors are 
arranged in descending order according to the F value. 

Significant effect on s/n: 
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Pre-etch line width: A, F 

Post-etch line width: None 
Significant effect on mean: 

Pre-etch line width: A, B, C 

Post-etch line width: B, G, C, A, H 
Significant factors identified by accumulation analysis: 

Post-etch window size: B, C, A, H, G 

Factors that have a significant effect on the s/n and the factors 
identified to be significant by the accumulation analysis are all control 
factors. Setting their levels equal to optimum levels niinimizes the 
process variability. Here the control factors are A, F, B, C, H, and G. 

To keep the process mean on target we use a signal factor. Ideally, 
the signal factor should have a significant effect on the mean, but 
should have no effect on the s/n. Then changing the level of the signal 
factor would affect only the mean. In practice, a small effect on the 
s/n may have to be tolerated. 

Among the factors (A, B, C, G, and H) that have a significant effect 
on the mean, factors A, B, and C are relatively strong control factors 
as measured by the F statistics for the accumulation analysis and the 
ANOVA for pre-etch line-width s/n. Also, these factors are relatively 
difficult to change during production. So A, B, and C are not suitable 
as signal factor. Between the remaining two factors, G and H, G has 
greater effect on the mean and also shows as a less significant factor in 
accumulation analysis. So exposure time was assigned to be the signal 
factor. 

The optimum levels for the control factors were selected as follows. 
The mask dimension (A) and the aperture (F) have a significant effect 
on the s/n for pre-etch line width. From Table V we see that the 
2.5-jum mask (level 2) has a higher s/n than the 2.0-jum mask. Hence 
2.5 jitm was chosen to be the optimum mask dimension. Also, aperture 
1 (level 1) has the highest s/n among the three apertures studied. 
However, because of the past experience, aperture 2 was chosen to be 
the preferred level. 

The accumulation analysis of the post-etch window-size data indi- 
cated that the viscosity, spin speed, mask dimension, developing time, 
and exposure have statistically significant effects on the frequency 
distribution. The optimum levels of these factors can be determined 
from Table XIII and Fig. 6 to be those that have the smallest traction 
of windows not open (category I) and the largest fraction of windows 
in the range 3.0 ± 0.25 [im (category IV). Because it is more critical to 
have all the windows open, when there was a conflict we took the 
requirement on category I to be the dominant requirement. The 
optimum levels are: 2.5-/zm mask dimension, viscosity 204, 4000-rpm 
spin speed, 60-second developing time, and normal exposure. 
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Table XV — Optimum factor levels 



Standard Optimum 
Label Factors Name Levels Levels 



A Mask Dimension (/on) 

B Viscosity 

C Spin Speed (rpm) 

D Bake Temperature (°C) 

E Bake Time (min) 

F Aperture 

G Exposure (PEP setting) 

H Developing Time (s) 

I Plasma Etch Time (min) 



2.0 


2.5 


204 


204 


3000 


4000 


105 


105 


30 


30 


2 


2 


Normal 


Normal 


45 


60 


13.2 


13.2 



Table XV shows side by side the optimum factor levels and the 
standard levels as of September 1980. Note that our experiment has 
indicated that the mask dimension be changed from 2.0 /xm to 2.5 jum, 
spin speed from 3000 rpm to 4000 rpm, and developing time from 45 
seconds to 60 seconds. The exposure time is to be adjusted to get the 
correct mean value of the line width and the window size. The levels 
of the other factors, which remain unchanged, have been confirmed to 
be optimum to start with. 

In deriving the optimum conditions we have conducted a highly 
fractionated factorial experiment and have considered only the main 
effects of the factors. The interactions between the factors have been 
ignored. If the interactions are strong compared to the main effects, 
then there is a possibility that the optimum conditions thus derived 
would not improve the process. So experiments have to be conducted 
to verify the optimum conditions. The verification was done in con- 
junction with the implementation, which is described next. 

X. IMPLEMENTATION AND THE BENEFITS OF THE OPTIMUM LEVELS 

We started to use the optimum process conditions given in Table 
XV in the Integrated Circuits Design Capability Laboratory in January 
1981. In the beginning the exposure was set at 90, which is the normal 
setting given in Table I. We observed that the final window at the end 
of step (xi) was much larger than the target size of 3.5 jum. Through 
successive experiments, we reduced the exposure time until the mean 
final window size came to about 3.5 jum. The corresponding exposure 
setting is 140. Since then the process has been run at these conditions. 
The benefits of running the process at these conditions are: 

(i) The pre-etch line width is routinely used as a process quality 
indicator. Before September 1980 the standard deviation of this indi- 
cator was 0.29 ixm on a base line chip (DSO chip). With the optimum 
process parameters, that standard deviation has come down to 0.14 
jum. This is a two-fold reduction in standard deviation, or a four-fold 
reduction in variance. This is evidenced by Fig. 7, which shows a 
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typical photograph of the PLA area of a BELLMAC-32 microprocessor 
chip fabricated by the new process. Note that windows in Fig. 7 are 
much more uniform in size compared to the windows in Fig. 2. Also, 
all windows are printed and open in Fig. 7. 

(ii) After the final step of window forming, i.e., after step (xi), the 
windows are visually examined on a routine basis. Analysis of the 
quality control data on the DSO chip, which has an area of approxi- 
mately 0.19 cm 2 , showed that prior to September 1980 about 0.12 
window per chip was either not open or not printed (i.e., approximately 
one incidence of window not open or not printed was found in eight 
chips). With the new process only 0.04 window per chip is not open or 
printed (i.e., approximately one incidence of window not open or 
printed is found in twenty-five chips). This is a three-fold reduction in 
defect density due to unopened windows. 

(Hi) Observing these improvements over several weeks, the process 
engineers gained a confidence in the stability and robustness of the 
new process parameters. So they eliminated a number of in-process 
checks. As a result the overall time spent by the wafers in window 
photolithography has been reduced by a factor of two. 

The optimum parameter levels were first used in the Integrated 
Circuit Device Capability Laboratory with only a few codes of ICs. 
Subsequently, these parameter levels were used with all codes of 3.5- 
jxm technology chips, including BELLMAC-A microcomputer and 
BELLMAC-32 microprocessor chips. The mask dimension change 
from 2.0 to 2.5 /xm is now a standard for 3.5-/rai CMOS technology. 

XI. DISCUSSION AND FUTURE WORK 

The off-line quality control method is an efficient method of im- 
proving the quality and the yield of a production process. The method 
has a great deal of similarity with the response surface method 9 and 
the evolutionary operations method, 10 which are commonly known in 
statistical literature in this country. Both the response surface and the 
evolutionary operations methods are used to maximize the yield of a 
production process and they both make use of the experimental design 
techniques. The main difference is that in the off-line quality control 
method the process variability that has a great impact on the product 
quality is the objective function. In the response surface and evolu- 
tionary operations methods, the process variability is generally not 
considered. Thus, intuitively, the optimum levels derived by using the 
off-line quality control method can be expected to be more robust, 
stable, and dependable. 

In the response surface method one typically uses a relatively large 
fraction of the factorial experiment. However, in off-line quality control 
usually a very small fraction is chosen. Another difference is that in 
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Fig. 7 — Typical results for contact windows in PLA area (January through March, 
1981) using parameters derived from this experiment. The contact windows are round 
shaped. 

the response surface method the objective function is considered to be 
a continuous function approximated by a low-order polynomial. In off- 
line quality control, we can simultaneously study both the continuous 
and discrete factors. 

Our application of the off-line quality control method to the window- 
cutting process in the Murray Hill 3.5-/mi CMOS technology, as seen 
from the earlier sections, has resulted in improved control of window 
size, lower incidence of unopened windows, and reduced time for 
window photolithography. Presently, we have undertaken to optimize 
two more steps in IC fabrication. Those steps are polysilicon patterning 
and aluminum patterning. Both these processes, like the window- 
cutting process, involve photolithography and are among the more 
critical processes of IC fabrication. We think that the method has a 
great potential and would like to see applications in various parts of 
Bell Laboratories and Western Electric Company. 
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APPENDIX A 

Computation of the Sum of Squares — Analysis of the Signal-to-Noise Ratio 
for Pre-Etch 

The computations of the sum of squares tabulated in Table Via are 
illustrated below. 

S m = Correction Factor 

18 

I* 

'- 1 ' (25.2202) 2 



= 35.3366 



18 18 

Sa = Sum of squares for factor A 
(9m Al ) 2 + (9m A f 



9 
(11.5711) 2 + (13.6491) 2 



— "3m 

- 35.3366 



9 
= 0.2399 (d.f. = 1) 

. (6m C] ) 2 + (6m c f + (6mc 3 ) 2 
Sc g S„ 



* Trademark of Bell Laboratories. 
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(8.1979) 2 + (8.1017) 2 + ( 8.9206) 2 
6 m 

= 0.0668 (d.f. = 2). 

Sums of squares for the factors E, F, G, and H were calculated 
similarly. The combined sum of squares due to B and D is given by 

Sbd = Sum of squares for the column BD 

(6m fllfll ) 2 + (6m Bl D 2 ) 2 + (ema^f 
6 bm 

= (S-2524) 2 + (8.6649) 2 + (8.3029) 2 _ 3& 3366 
6 

= 0.0169 (d.f. = 2). 

The total sum of squares is 

S T = £ tj 2 - S m - 0.8423. (d.f. = 17). 

The error sum of squares is calculated by subtraction. 

S e = S T - (S A + Sbd + S c + Se + Sf + S g + S h ) 

- 0.1522 (d.f. = 4). 

Here we do not compute the sum of squares due to factor I (etch 
time), because it has no influence on the pre-etch line width. 

APPENDIX B 

Computation of the Percent Contribution — Analysis of the Signal-to-Noise 
Ratio for Pre-Etch Line Width 

The computation of the percent contribution is explained below. 
The contribution of factor A to the total sum of squares 

= S A - (d.f. of A) (error mean square). 

Hence, the percent contribution for factor A 

Sa — (d.f. of A) (error mean square) 

= X 100 

total sum of squares 

,0.2399 -0.0281 

0.8423 

The percent contributions of E and F are determined similarly. Now 
consider, the contribution of error to the total sum of squares: 

= S e + (total d.f. for factors) (error mean square). 
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Hence, the percent contribution for error 

S e + (total d.f. for factors) (error mean square) 

total sum of squares 
0.3010 + 5 X 0.0251 



X 100 



0.8423 



X 100 = 50.6%. 



APPENDIX C 

Separation of Sbd into S B and S D — Analysis of the Mean Pre-Etch Line Width 

The sum of squares, Sbd, can be decomposed in the following two 
ways 11 to obtain the contributions of the factors B and D: 

Sbd = S'b(D) + Sb 

and 

Sbd = Sb(B) + S'b. 

Here S'b(d) is the sum of squares due to B, assuming D has no effect; 
Sb is the sum of squares due to D after eliminating the effect of B. 
The terms Sb(B) and S'b are interpreted similarly. We have 

S B{D) <i« + i» + 2 »)x6 0.903 (d.f. - 1) 

Sb = Sbd ~ S' B (D) = 0.047 (d.f. = 1). 
Similarly, 

o, _ (QmB lDl + Qm Bi p x - 12m Bl D 2 ) 2 _ ( . f 

S D{B) (1 . + l2 + 2 2 )x6 °- 116 W- - « 

Si = Sbz> - Si,( B) = 0.834 (rf./. = 1). 

For testing the significance of the factors B and D we use S'b and 
S' D , respectively. Note that S' B and Sb do not add up to Sbd, which is 
to be expected because the design is not orthogonal with respect to the 
factors B and D. 

APPENDIX D 

Computation of the Sum of Squares for Accumulation Analysis — Analysis of 
Post-Etch Window Size 

The weights for the cumulative categories (I), (II), (III), and (IV) 
are given below. The frequencies of the bottom line of Table XIII are 
used in computing these weights. Therefore: 

W<i) =-— ^r = 4.008 

86 180 - 86 

180 X 180 
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180 
W m = 



Wan) = 



102 x (180 - 


% - 4.072 
•102) 


180 2 


— a 713 


125 X (180 - 


• 125) 


180 2 


= 10.125. 



WilV) 160 X (180 - 160) 

Computation of the sum of squares tabulated in Table XIV are 

illustrated below: 

/51 2 + 35 2 86 2 \ „, /58 2 + 44 2 102 2 
Sa = W m X t^ + Won X 



90 180; V 9° 180 

/67 2 + 58 2 125 2 \ _, /88 2 +72 2 160 2 

+ W m) x — -— r + WiivjX 



90 180 ; \ 90 180 

= 26.64, 

and 

/47 2 + 22 2 + 17 2 86 2 \ 
S C =W (I) X^ — J 

/52 2 + 29 2 + 21 2 102 2 
+ VVl " lX| 6^ MO 

58 2 + 39 2 + 28 2 125 2 



+ Warn x 
+ Wav) X 



60 180 

60 2 + 55 2 + 45 2 160 2 



60 180 

= 125.52. 
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